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Abstract 

Some of the consequences for inflationary cosmology of a scale dependence 
(running) in the tilt of the scalar perturbation spectrum are considered. In 
the limit where the running is itself approximately scale-invariant, a relation- 
ship is found between the scalar and tensor perturbation amplitudes, the scalar 
spectral index and its running. This relationship is independent of the func- 
tional form of the inflaton potential. More general settings, including that 
of braneworld cosmological models, are also considered. It is found that for 
the Randall-Sundrum single braneworld scenario, the corresponding relation 
between the observables takes precisely the same form as that arising in the 
standard cosmology. Some implications of the observations failing to satisfy 
such a relationship are discussed. 
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1 Introduction 



The Wilkinson Microwave Anisotropy Probe (WMAP) has measured the power spec- 
trum of the Cosmic Microwave Background (CMB) for multipoles up to I ~ 800 with 
unprecedented accuracy p. The best-fit model to the WMAP data alone is consis- 
tent with a spatially flat universe, with near scale-invariant, adiabatic and Gaussian 
distributed primordial density (scalar) perturbations [3 El El, as predicted by the 
simplest models of inflation [5 . (For a review, see, e.g., Ref. [0]). 

Combinations of CMB results with other astrophysical observations have led to 
strong constraints on the standard cosmological parameters such as the Hubble pa- 
rameter, baryon density and age of the universe j2] . In order to differentiate between 
the numerous inflationary models, however, it is necessary to constrain the power 
spectrum of the primordial fluctuations. Assuming that such a spectrum varies in a 
suitable fashion, the standard approach is to expand its logarithm as a Taylor series 
in Ink, about a given scale, k^. 

In A s (k) = In A s (ko) + (n s - 1) In A + J-a 5 In 2 £ + . . . , (1) 

k 2 k 

where k is the comoving wavenumber, ns is the spectral index (tilt) of the spectrum 
and the second-order term, as = (dns/dlnk)^ , represents the 'running' of the spec- 
tral index [7 . The 'power-law' approximation is equivalent to truncating the spectrum 
to first order, i.e., specifying as = 0. At this level of approximation, the best-fit to the 
WMAP data is ns = 0.99 ± 0.04 j2J. On the other hand, there is some evidence that 
the power-law approximation may be inadequate when data sets spanning a much 
wider range of scales are combined. Specifically, Peiris et al. |3] include CMB data 
from the CBI |H] and ACBAR (covering the range of multipoles 800 < I < 2000 
complementary to WMAP), together with the two degree field (2dF) galaxy redshift 
survey JU] and Lyman-a forest data at wavenumbers above k ~ 0.1 Mpc [TT]. There 
is marginal 1.3a support for a non-zero running, as = —O-O^-oml 13- However, the 
validity of employing the Lyman-a forest data has been questioned [12] . Bridle et al. 
[T3] include CMB data from the VSA [T4J but not from the Lyman-a forest, and find 
the marginalised la result as = —0.04 ± 0.03, in agreement with the WMAP collab- 
oration, although they conclude that evidence for a non-zero running is dependent 
on the surprisingly low values of the quadrupole and octopole moments in the CMB 
power spectrum. In particular, as = is consistent when the I < 5 multipoles are 
excluded [THj. Other authors who include only CMB and 2dF data also find that a 
scale-invariant tilt is consistent with the observations fSl EE El ■ 

Although there still remain some open questions regarding the interpretation of 
the data, the recent developments outlined above provide strong motivation for con- 
sidering what one might expect to learn about inflationary cosmology if a running 
in the spectral index is detected HH1 [HI [201 EH 122] • This is especially true given 
the anticipated improvement in the quality of data from future satellite experiments 
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such as Planck. In general inflationary settings one would expect the running itself 
to be scale-dependent. However, given the current absence of observational evidence 
for such a variation and, furthermore, that a varying running may be approximated 
as a piecewise constant over a small enough range of scales, we consider inflationary 
models where the running of the tilt is scale-independent and non-zero. It is found 
in Section 2 that the amplitude, tilt and running of the scalar spectrum are related 
in a non-trivial fashion to the amplitude of the gravitational wave spectrum that is 
also generated during inflation. A similar relationship is found in Section 3 for a class 
of braneworld inflationary cosmologies. We conclude with a discussion in Section 4. 



2 Running of the Spectral Index 

In general, the power spectrum of the scalar perturbations is closely related to the 
functional form of the inflaton potential 1 , V(<p): 

where n 2 = 87rmp 2 , mp is the Planck mass and a prime denotes d/d(j). The rela- 
tionship between the inflaton field and comoving wavenumber follows from the scalar 
field equations of motion and is given by 

" V ' " (3) 



(4) 



d\nk 3H 2 d(j) 
in the slow-roll limit. By defining the 'slow-roll' parameters as [25J: 

1 V' 2 

e= 2^ 2 V I 
V" 

V'V" 

the spectral index and its running may be expressed directly in terms of the potential 
and its derivatives [7| l2l)j: 

n s - 1 = -6e + 2r] (7) 
a s = I6er] - 24e 2 - 2f . (8) 

Thus, the running of the spectral index depends on the third derivative of the poten- 
tial. Eq. is truncated at order £, such that quadratic corrections in e and r\ are 

^^We employ the normalization conventions of Ref. |24| in this paper. 
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assumed to be negligible. This requires that |£| <C max (e, \r]\) and is equivalent to 
assuming that |ng — 1| <C 1 and \ag\ ~ (ns — l) 2 or less. As emphasized in Ref. |19j . 
slow-roll predicts the former condition but not necessarily the latter. 

The slow-roll parameter (@J) is also related to the spectrum of tensor (gravitational 
wave) perturbations that are generated quantum mechanically during inflation |23j . 
The relationship is expressed through the consistency equation (for a review, see, e.g., 
Ref. j21): 

A 2 1 

= ~2 nT ' Ut = _2e ' ^ 

s 

where {At, ut} represent the amplitude and spectral index of the tensor perturbations 
respectively. 

In view of the discussion given in the previous Section, we consider the case where 
as is assumed to be constant. Our aim is to derive an expression relating observ- 
able parameters in the presence of a non-zero running. This requires the integration 
of Eq. (JHJ) with respect to the inflaton field. This equation may be viewed as a 
third-order, non-linear differential equation. Its first integral therefore relates the 
inflaton potential to its first two derivatives, or equivalently, the two slow-roll pa- 
rameters (jlj) and (j5J). Consequently, substitution of Eqs. (J7J) and (J0J) into such an 
expression then results in a constraint equation that relates the observable parameters 
{A s ,n s ,as,A T }. 

To proceed, we introduce a new variable 

V 

y = y (10) 

representing the logarithmic derivative of the potential. The third-order equation (JSJ) 
then reduces to the non-linear, second-order equation 

4 

yy -y y = — ^— ■ (n) 

Defining z = y', such that y" = zdz/dy, then reduces Eq. (II lj) to a first-order 
equation of the form 

dz 2 K 4 a s no s 

yz—-yz= — . 12 

dy 2 

Eq. (|12j) may be further simplified by defining the variable 

u = z- l -y 2 (13) 
and it follows after substitution of Eq. (|T3^) into Eq. (|T2^) that 



dy = 2_ 

du K 4 a>s 



1 3 

uy + -y 



(14) 
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where we now view y and u as the dependent and independent variables, respectively. 
Eq. (|14j) may then be rewritten in a separable form by introducing the variable 



w = y(u) exp 



u 



(15) 



and substitution of Eq. (JT5j) into Eq. (fT^j) implies that 



1 dw 

w 3 du 

Eq. (JTfij) admits the first integral 



K 4 as 



exp 



2u 2 



(16) 



exp 



2m 2 



7T 



2n A a s 



1/2 



erf 



(17) 



where erf(x) = 27r~ 1 / 2 Jq ds exp[— s 2 ] represents the error function and c is an in- 
tegration constant. The error function is a monotonically increasing function, such 
that erf(O) = and erf(oo) = 1, and has a first derivative given by d [erf(x)] jdx = 



(2/ x /vr)exp( 



-x 



In the case where a,s < 0, Eq. (fTTjl may be expressed in terms of 



the imaginary error function, erfi(a;) = — ievf(ix). 

Finally, the pair {u, y} may be directly related to observable parameters. Compar- 
ison of Eqs. (@J) and (fTUJl implies that y 2 = 2n 2 e and it then follows from definition (fT3j) 



that the variable u is directly related to the scalar spectral index, u 
Thus, substituting Eq. (J§J) into Eq. (fTTjl implies that 



(n 5 -l)/2. 



A 2 

A, 



2 exp 



2a s 



O^x 1/2 
Z7T \ ' 

a s 



erf 



n s 



V2^s 



where c is an undetermined constant that is, in principle, measurable. Eq. ()18|) rep- 
resents an observable signature of inflationary models that generate a scale-invariant 
running of the spectral index. 



3 Running and Braneworld Cosmology 

It is also of interest to investigate whether the above type of observable signature holds 
in other more general inflationary scenarios. In recent years, considerable interest 
has focused on the possibility that our observable four-dimensional universe may be 
viewed as a domain wall or 'brane' that is embedded in a higher-dimensional 'bulk' 
space [23l2Ell2n]- According to these scenarios, the standard model gauge interactions 
are confined to the brane, but gravity may propagate in the bulk [28J. The motion of 
the brane through the static bulk space is interpreted by an observer confined to the 
brane as cosmic expansion or contraction |3D] . 



4 



In this Section, we consider the Randall-Sundrum type II (RSII) braneworld sce- 
nario, where a single brane is embedded in five-dimensional anti-de Sitter (AdS) 
space [21]. In this case, the effective Friedmann equation on the brane is derived 
from the Israel junction conditions that relate the extrinsic curvature of the induced 
metric on the brane to the energy-momentum of the matter fields that are confined 
to the brane [31] . The form of the Friedmann equation is modified from that of stan- 
dard cosmology based on Einstein gravity and acquires a quadratic dependence on 
the energy density, p [32] : 

(19) 



1 + £ 



where A is the brane tension. 

Such a modification becomes important at high energies and has significant impli- 
cations for inflation [331 EH]- m particular, in the limit where p ^> A, the amplitudes 
of the scalar and tensor perturbations are enhanced 



3 ~ BOOttW 2 ' t ~ 200^' ( ' 

However, despite these corrections, the consistency equation relating the two spectra 
is identical to that of the standard scenario, Eq. 36 j. Such a degeneracy in the 
consistency equation also arises in more general brane cosmologies |37j . 

A natural question to address, therefore, is whether such a degeneracy may be 
lifted by allowing for a running of the spectral index. In view of the modifications to 
the Friedmann equation that typically arise in brane cosmology, we consider a model 
described by a Friedmann equation of the form 

H 2 = (21) 

where p is the energy density of the matter, q is an arbitrary, positive constant and 
k 2 is an arbitrary constant. Eq. (j21[) may be viewed as a limiting case of a more 
generalized Friedmann equation that is relevant in the high energy regime of early 
universe dynamics. For example, the case q = 2 corresponds to the RSII scenario 
when the energy density dominates the brane tension and we specify k 2 = k 2 /2X. A 
further case of interest is given by q = 2/3. A Friedmann equation of this form arises 
in the extended version of the RSII scenario when a Gauss-Bonnet combination of 
curvature invariants is included in the five-dimensional bulk action [35] . 

We further assume that the universe is dominated by a single, self-interacting 
inflaton field. Conservation of energy-momentum of this field then implies that 

<j> + 3H0 + V = 0, (22) 

where a dot denotes differentiation with respect to time. We define the generalized 
slow-roll parameters as e g = —H/H 2 , rj g = V"/(3H 2 ) and £ g = V'V" /(3H 2 ) 2 , 



5 



respectively, and in the slow-roll limit, <p 2 <^V and \(f>\ <C H\<p\, these reduce to 



e 9 - 2R2 vq+ x y Z6 ) 

V" 

% = SPPI (24 > 

V'V" 

^ = ~iyV ( 25 ) 

Conservation of energy-momentum implies that the curvature perturbation on 
uniform density hypersurfaces is conserved on super-Hubble radius scales. This fol- 
lows as a direct consequence of energy-momentum conservation of the inflaton and is 
independent of the gravitational physics (HHl- It can then be shown that the amplitude 
of the scalar perturbation spectrum is given by Ag oc H /(f) 2 in the slow-roll limit 
|3*9"] . The value of the scalar field is related to the comoving wavenumber through Eq. 
(JHJ)- Substituting the field equations into the expression for the amplitude and differ- 
entiating with respect to comoving wavenumber then implies that the scalar spectral 
index is given by 

n s -l = -6e 9 + 2 % (26) 
and the running of the tilt is given by 

as = 16e^ - ^±^-e 2 g - 2i r (27) 

As in the previous Section, our aim is to integrate Eq. (J2*7)l under the assumption 
that the running of the spectral index is constant. It proves convenient to define a 
new scalar field, tp: 

d _ V (i- q )/2d 
dip d(p 
and this implies that Eq. (|2"Tj) takes the form 

K*a s = A{q + 1)^— - 2{2q + 1)VLL - 2^^, (29) 

where a star denotes d/dip. Defining the new variable Y = V*/V then simplifies Eq. 
flU to 

~4 

YY** + (1 - 2q)Y 2 Y* = -1^1. ( 30 ) 

Eq. (fHU|) can be integrated in a similar way to that employed in Section 2 and we 
therefore omit the details. Eq. (|HU|) reduces to the separable equation 



1 dW (l-2q 



W 3 dU V k 4 a s 



exp 



2U 2 

k 4 as 



(31) 
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where 



U = Y* + 



'1-2^ 


Y 2 

2 ' 


, 2 
exp 


) 

' U 


k 4 as 



(32) 
(33) 



and solving Eq. (}3*Tj) then implies that 



1 

F2 



exp 



" 2U 2 ~ 













7T 



2£ 4 a s 



(1 - 2g)erf 



U 



c, 



(34) 



where C is an integration constant. 

Comparison of Eqs. ()26|) and ()32|) implies that U = k 2 (ns — l)/2. Moreover, 
substituting Eq. (J2*%|) into the definition of the variable Y implies that Y 2 = (2k 2 / q)e g . 
It follows, therefore, that Eq. (JHlj) may be expressed in the form 



1 



exp 



{ns 



2«c 



,a s J \ q J 



C. 



(35) 



where C is a dimensionless constant. 

In braneworld inflationary scenarios, the calculation of the tensor perturbation 
spectrum is more involved than that of the scalar perturbations because the gravita- 
tional waves extend into the bulk dimensions [S3]. Consequently, one must consider 
the tensor perturbations for each specific model. For the RSII scenario, where q = 2, 
substituting Eq. (jZU) into Eq. fl2D|) implies that A 2 T /A 2 S = 3e g /2. Remarkably, 
therefore, we conclude that when Eq. (|33|) is expressed in terms of the observables 
{As, ns, as, At}, it reduces to precisely the same form as the corresponding relation- 
ship for the standard inflationary cosmology, Eq. (JTHJ). 



4 Discussion 

The inflationary scenario has received a great deal of observational support from re- 
cent CMB satellite observations From the theoretical perspective, an important 
problem to address is the origin of the inflaton field within a fundamental underlying 
theory and, more specifically, the nature of the inflaton potential that drove the ac- 
celerated expansion of the very early universe. In the case where the running of the 
spectral index vanishes, it is well known that the form of the potential leading to a 
constant spectral index is not unique j2HH0]. Indeed, the origin of this degeneracy 
may be understood from a mathematical point of view by expressing the potential in 
terms of the derivative V = dW ($) / d(f> and rewriting Eq. (J7J) in the form 



W"' 




(36) 
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This novel way of expressing the constraint on the potential is illuminating because 
the left hand side of Eq. (J36j) is the Schwarzian derivative of the function W (0) 
[4~T] . This is the unique elementary function of the derivatives that is invariant under 
the homographic transformation that corresponds to the group of fractional linear 
transformations: 

cW + d K ' 

where {a, b, c, d} are arbitrary constants satisfying ad — be = 1. Thus, given a partic- 
ular solution to Eq. (j3T?j) (such as an exponential potential), more general solutions 
and corresponding potentials may be generated by applying the transformation (J37j) . 
Further observational input, most notably from the gravitational wave background, 
is required to lift the degeneracy 

In this paper we have considered the more general class of inflationary models 
where the running of the scalar spectral index is itself non-zero and independent of 
scale. In general, it is not possible to determine the analytic form of the potentials 
that generate such a spectrum. On the other hand, their asymptotic limit may be 
deduced by noting that the first integral of the second-order equation (JTTj) may be 
written in the form 



</' = ^f(<j>)-rfa s hiy, (38) 
where the function /(</>) itself satisfies the non-linear equation, /' = 2yy' 2 , i.e., 

/' = 2y(f -K A a s In y) . (39) 



The pair of equations (J58)l and (|3"9*|) represent a plane autonomous system with a 
single equilibrium point located at / = k as In y and it follows from the definition of 
the function / that this point represents the asymptotic form of the general solution 
to Eq. (fTT|) in the limit where \y"\ <C y\y'\, i.e., where the first term on the left hand 
side of Eq. ffTTj) is negligible. In this limit, Eq. (fTTj) may be integrated to yield the 
form of the potential: 



V = Vq exp 



(40) 



where Vq is an arbitrary positive constant and the sign of the exponent corresponds 
to the sign of the running. It is worth remarking that potentials of this specific 
asymptotic form also arise within the context of supergravity models [42j . 

Eq. (fT8|) directly relates the four observable parameters {As, ns, as, At}- An 
important feature of this relation is that it is independent of the specific functional 
form of the inflaton potential. In this sense, therefore, it represents an observable 
signature of the class of inflationary models where the 'running of the running', (3s = 
dots I d In k, is negligible. However, the expressions for the tilt and running, Eqs. ((7j) 
and ©, were derived under the standard assumption that |£| -C max(e, |^|). Although 
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this is consistent, it is not required by the slow-roll approximation ^Hj arid failure 
to satisfy Eq. (jTHj) may therefore indicate that such an assumption would need to be 
relaxed. 

The occurrence of the arbitrary constant implies that to proceed observationally 
the parameters {As, ns, as, At} must be measured over at least two separate scales. 
One measurement is required to determine the numerical value of the integration 
constant and the second to determine whether Eq. (jTHjl is indeed satisfied. The 
advantage of Eq. (jTHj) over the consistency equation Q is that it relates the scalar 
and tensor perturbation amplitudes directly to the scalar spectral index and its run- 
ning. Consequently, it does not require the tilt of the tensor spectrum to be directly 
measured. 

On the other hand, we have assumed that the running is effectively constant and 
this can only be verified observationally to within some error. When discussing obser- 
vational constraints, the primordial power spectrum may be viewed as an unknown 
function and the fitting procedure effectively truncates the Taylor expansion (JTJ) at 
a finite order. This is equivalent to setting all corresponding higher-order, slow-roll 
parameters to zero. Establishing whether a constant running is a good fit to the data 
requires the introduction of the next-order term, (3s = dots/ dink, as an additional 
parameter in the analysis. Self-consistency of the assumptions made above requires 
that f3s ~ \ns — 1| 3 or smaller. The assumption that the running is constant would 
then be supported if it turned out that (3s has only a moderate influence on the likeli- 
hood distributions for the other observable parameters and is itself entirely consistent 
with zero within the observed errors. 

Nevertheless, Eq. (jTHj) may prove important even if a high running of the running 
is reported. Leach and Liddle have argued that appropriate conditions should be 
satisfied if the inclusion of a higher-order parameter of the power spectrum is to be 
justified jT7j . In effect, the criterion is that of convergence in the Taylor expansion 
(JJJ. In the present context, the inclusion of the running of the running could only be 
justified if the third-order term in Eq. (JTJ) is significantly smaller than the second- 
order term and, quantitatively, this requires 



If condition (j41|) is violated when the detection of (3s has only a low significance, it 
could be argued that the determination of this parameter may be unreliable and that 
it is therefore not appropriate to include it in the analysis [TZj. Indeed, Eq. (jTHjl proves 
important in this case because it may be employed to yield crucial information about 
the magnitude of the second derivative of the tilt. Failure of Eq. (jTHj) to satisfy the 
data (under the assumption of a constant running) could be interpreted as evidence 
that the second derivative of the tilt is important without the need for allowing this 
parameter to be a priori non-zero. 

We have also considered inflationary models that generate a scale-invariant run- 




(41) 
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ning of the spectral index in the Randall-Sundrum type II braneworld scenario. Sur- 
prisingly, we found that in the high energy limit, the relation (J35|) . when expressed 
in terms of the observables {As, ns, as, At}, takes precisely the same form as the 
corresponding relationship for the standard inflationary cosmology, Eq. (fTBj) . This 
provides further evidence of the degeneracy that exists between the primordial pertur- 
bations that are generated in the two scenarios even though the gravitational physics 
is manifestly different in the two cases |3T?1 13*7] . It should be borne in mind, however, 
that in these calculations the effects of the bulk space on the evolution of the density 
perturbations has been neglected. This is consistent at linear order when consid- 
ering scalar perturbations of a homogeneous background [53] . More generally, the 
backreaction perturbs the bulk space away from conformal invariance and generates 
a non-trivial Weyl curvature in the bulk WQ . This plays the role of a non-local 
energy-momentum source when projected down to four-dimensions and thus alters 
the background dynamics [I3JEI1E2I- The failure of the relation (JTHJ) to be satisfied 
in this case could therefore indicate the possible importance of these bulk effects. 
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